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We consider a disordered superconducting film in tlie magnetic field close to the upper critical 
field. Assuming that the phase transition to the glassy superconducting state exists and it is of the 
second order we show that it is exactly described by a zero-dimensional replica model having a form 
of one describing the phase transition in the SK (Sherrington-Kirkpatrick) spin-glass model. The 
dependence of the magnetic moment on the external magnetic field is shown to be more smooth 
than the Abrikosov mean field result for the clean case. While the total magnetic moment of the 
film close to the phase transition is diamagnetic, the part of it emerging due to the phase transition 
is paramagnetic, this shows tendency to the paramagnetic Meissner effect. 
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A magnetic field penetrates type-II superconductors 
forming the Abrikosov lattice of magnetic verticesEJ. Ob- 
viously, formation of a purely periodical structure is pos- 
sible only in homogeneous materials. The influence of 
disorder on the structure of the Abrikosov lattice was 
first considered in Reffl assuming that the random force 
acting on the lattice is uncorrelated with the lattice dis- 
placements and it was shown that the long range order 
is destroyed by inliomogeneity. Later this problem was 
considered in Ref.cl using more complicated ( but still not 
rigorous) methods where the result of absence of the long 
range order was confirmed and the correlation functions 
were analyzed more carefully. It, is expected theoretically 
and it was seen experimentallyQ that this system should 
have glassy properties. Therefore a systematic treatment 
of this problem is very complicated and is absent at the 
present time. 

We consider a disordered superconducting film in the 
external magnetic field close to the upper critical field 
Hc2 and show that the phase transition to the glassy 
superconducting state is exactly described by a zero di- 
mensional replica model having a form of one describing 
the phase transition in the infinite-range SK spin-glass 
model. To obtain this result we assume that the phase 
transition exists and it is of the second order. The depen- 
dence of the magnetic moment M on the external mag- 
netic field (or on temperature) following from this theory 
is smooth, only the first derivative dM/dH (dM/dT) has 
a kink in contrast with the Abrikosov mean field result 
for the clean case. The part of the magnetic moment 
emerging due to the transition to the glassy supercon- 
ducting state is paramagnetic, i.e. the system exhibits 
tendency to the paramagnetic Meissner effect. 

Our starting point is the Ginzburg-Landau free energy 
density 
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where D = V -I- ie*A/c, e* = 2|e|, and Ho is an ex- 
ternal magnetic field which is perpendicular to the film. 



The order parameter field (p in (^ is supposed to be 
normalized such that a = T — T^, where T^ is the 
mean field critical temperature. The width of the film 
d is supposed to be smaller than the coherence length. 
The potential u{r) models the disorder in the system 
and it is assumed to be Gaussian with the correlator 
(u(ri)u(r2)) = y/u S{ri — r2). This term emerges because 
in the random system the coefficient a = T — Tc becomes 
a random quantity. In fact all the coefficients in (m) con- 
tain randomness, but apparently, close to the phase tran- 
sition the most relevant effect comes from randomness of 
the coefficient a which we consider. 

Using the replica trick we can average (0) over u{r) 
getting the effective functional 
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where a, P are the replica indexes taking values from 1 to 
n. In the lowest Landau level approximation the fluctu- 
ations of the magnetic fleld can be integratiid out giving 
the renormalization of the interaction termtl 
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H^/Stt. We 



where v = b— {e*/2mc)^/2n and T' = T 
shall take the gage A^, — —Hy and use the Landau func- 
tions to diagonalize the quadratic part of the functional 
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where I = ^c/e*H , L^ is the size of the film in the 
X— direction, and p is the momentum. In this basis the 
functional (S) becomes 
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At the glass transition this term becomes zero. This gives 
the value of g at the phase transition 
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where 5 = {e*H/mc + a)/Tc and the kernel K is 
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Note that the Lagrangian ds) is translationally invariant 
in the momentum p. In the clean case this invariance 
is broken when the Abrikosov lattice is formed. But in 
the disordered case which we consider we expect that 
all averaged quantities are translationally invariant in p. 
Therefore the averaged Green function 
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should not depend on p. The functional (||) as usually can 
be presented as an expansion in the total Green function 
Ga^pip)- Deriving this effective functional one can per- 
form all momentum integrations in each diagram because 
the momentum dependence of the Green function is triv- 
ial, and the resulting free energy functional contains Ga,p 
depending only on the replica indexes. Keeping the terms 
up to the forth order in G we get 
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where the dimensionless coupling constants are k — 
v/inTcPd, 6 = u/inTaPd, and N = e*HoLxLy/2'Kc is 
the number of quantum states at the lowest Landau level. 
Let us first analyze the model (pi) and then we will dis- 
cuss how the results are affected by the high order terms 
neglected in (||). Variation of (g) with respect to Ga.p 
gives the equation defining Ga.p- In the normal phase the 
solution should be replica symmetric Ga^j3 = g5a,i3 be- 
cause the model (0) is symmetric under the global C/(l) 
transformations within each replica. Variation of (^ with 
respect to g gives 

- - + 5 + 2(k - 6l)g + (4k6' -k^- 29^)g^ = (9) 

In the glassy superconducting state we expect a non- 
replica- symmetric solution for Ga,f3 to appear. The glass 
transition corresponds to the instability of the replica- 
symmetric solution. Presenting Ga,i3 as a sum of diag- 
onal and off-diagonal parts Ga^p — gSa^p + Qa,f} and 
expanding (0) to the second orderer in Q we get 
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We see that within the approximation which we used the 
second order phase transition is possible only when the 
disorder is strong enough 0/ k > 4/7. The value of 6 at 
the phase transition is determined by (0) 
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We will be interested in the region close to the glass tran- 
sition where we can simplify (0) expanding in Q. Taking 
g = gc+g' , 5 = 5c+5' , expanding in Q,g' and eliminating 
g' we get the following effective free energy functional 
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We do not write the expressions for the coefficients a2, as 
because they do not affect the results in the leading 
orderS. The model (^3|) is standard in the spin-glass the- 
ory, it describes the phase tcansition in the SK spin-glass 
model (see for example ReLQ). The solution is known to 
have the form of the Parisi matrix which is parameterized 
by the function q{x), < a: < 1, after taking the limit 
n -^ 0. For the model (O) the function q{x) is 
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where x„ 



-^5' and q^ = -'-. 
Now let us analyze the validity of the obtained results: 
According to ( |ll| ) the Green function Ga^p is of the or- 
der of \/\9. Therefore the high order corrections to (j^) 
are of the same order or less than the terms that were 
kept. It means that the values of the coefficients (|lj, 
[l5| , [l6|) are correct only by the order of magnitude. Also 
Eq. (|1 1[) can be used only as an estimation, therefore the 
existence of the second order phase transition in fact is 
an assumption. Nevertheless, the form of the functional 
( |l3| ) is correct because it contains all terms up to the 
forth order allowed by the replica symmetry of the origi- 
nal model and the expansion in Q is always possible close 



to the phase transition. Therefore the phase transition 
is effectively described by a zero-dimensional model (not 
including the replica indexes). While this is typical for 
the infinite range spin-glass models it is a very surprising 
result for a real physical system. 

The magnetic moment of the system can be written as 
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where Mq is the part of the magnetic moment emerging 
due to formation of the glass state 
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and Mfi is the part which can not be found from the 
effective functional (|l^) but it should not have any pe- 
culiarities at the phase transition. In the normal state 
Mq — and in the glass state from (O) wc have 
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We see that the dependence of the magnetic moment 
on the external magnetic field is smooth only the sus- 
ceptibility dependence dM/dHo has a kink at the phase 
transition. Since the temperature enters the functional 
(o) in the same way as the magnetic field, analogously 
we conclude that the entropy is a smooth function of 
the temperature and the specific heat dependence on the 
temperature has a kink. For a comparison let us note 
that in the pure case the Abrikosov mean field theory 
predicts more singular behavior: the magnetic moment 
and specific heat dependences should have kinks at the 
phase transition. 

Within the lowest Landau level approximation the to- 
tal magnetic moment (|l9|) is always negative, this is clear 
from (^ because {4'a.p4'a,p) > 0. But the part of the mag- 
netic moment emerging due to the transition to the glassy 
superconducting state is positive according to (|2l|). It 
means that the system has tendency to exhibit the para- 
magnetic Meissner effect. 

This theory can be applied for experiments on Nb su- 
perconducting films in magnetic fields close to Hc2- To 
the best of our knowledge the region close to the phase 
transition was not yet investigated. Our theory predicts 
classical critical exponents for dependences of magne- 
tization and specific heat on temperature or magnetic 
field. The dependence of the magnetization M on the 
temperature or magnetic field is smooth, only the first 
derivatives dM/dT, dM/dHo should have kinks. Since 
our theory predicts that the magnetic moment emerging 
due to the phase transition is paramagnetic, experimen- 
tally one would expect that decreasing the temperature, 
the magnetization first decreases due to the effect of the 
diamagnetic fluctuations, but then increases bellow the 
phase transition. 

It is interesting to^ mention the experiment on the su- 
perconducting diskgj where the phase transition region 



was investigated in details. The results of this experi- 
ment qualitatively agree with our theory ,but the thick- 
ness of the disks that were used is larger then the coher- 
ence length, therefore our theory can hardly be directly 
applied. 

I am grateful to Elihu Abrahams, Thierry Giamarchi, 
Lev loffe, Gabriel Kotliar and A.I. Larkin for very useful 
discussions. 
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